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Abstract 

We present a version of acoustic black holes by using the principle of the Josephson effect. We 
find that in the case two superconductors A and B are separated by an insulating barrier, an acoustic 
black hole may be created in the middle region between the two superconductors. We discuss in 
detail how to describe an acoustic black hole in the Josephson junction and write the metric in the 
langauge of the superconducting electronics. Our final results infer that for big enough tunneling 
current and thickness of the junction, experimental verification of the Hawking temperature could 
be possible. 
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1 Introduction 

Analog models of gravity have recently received great interests since these models may provide possible 
experimental verifications of the evaporation of black holes. Unruh was the first to propose the idea of 
using hydrodynamic flows as analogues to mimic some properties of black hole physics HI (for reviews 
see and references there in). Any moving fluid with speed exceeding the local sound velocity 
through a spherical surface could in principle form an acoustic black hole. For acoustic black holes, it 
is sound waves instead of light waves that cannot escape from the horizon where the horizon locates 
on the boundary between subsonic and supersonic flow regions. In particular, superfluid helium II 
1131 , atomic Bose-Einstein condensates |4]-(7]|, one-dimensional Fermi-degenerate noninteracting gas (8|| 
were proposed to create an acoustic black hole geometry in the laboratory. 

Because the Hawking temperature depends on the gradients of the flow speed at the horizon, de- 
tecting thermal phonons radiating from the horizons is very difficult. In fact the Hawking temperature 
calculated from models in Bose-Einstein condensates so far is very low (~ nano Kelvin). Up to now, 
only few experimenters have claimed that acoustic black holes are able to be observed. The very recent 
experimental realization of acoustic black hole reported was conducted in a Bose-Einstein conden- 
sate [9|. Considering that the Bose-Einstein condensate systems have a very strict requirement on the 
environment temperature (for example 170-nano Kelvin for the gas of Rubidium atoms), the authors 
in [[T0]| proposed that acoustic black holes may be realized in superconducting materials which have 
much higher critical temperatures (~ 100 Kelvin) and a relativistic version of acoustic black holes was 
presented there (see also ifTTT - fTol for further reading). It was observed that an acoustic black hole may 
form near the spiral vortex core in a type II superconductor iTTOl . But the experimental detecting of 
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such behavior in superconductors is so far very difficult, because the superconducting phase could be 
destroyed when the speed of the current carriers exceeds the Fermi-velocity. The purpose of this paper 
is to investigate the possibility of creating acoustic black holes in the supercurrent tunneling, in partic- 
ular the Josephson junctions, because there are some very sensitive detectors based on the principle of 
the Josephson effect, for example, the superconducting quantum interference device (SQUID). 



2 Brief reviews on the Josephson effect 

In the theory of superconductivity, the Ginzburg-Landau equation is a phenomenologically based the- 
ory, which has been proven to be very successful partly because it can describe the mixed states in type 
II superconductors. It was later shown by Gorkov that this equation can be derived from the Bardeen- 
Cooper-Schrieffer (BCS) theory in the case of thermodynamic equilibrium and for temperatures close 
to the transition temperature IfTTl . The generalization of the Ginzburg-Landau equation, which includes 
the time dependent term, allows the superconductor to relax in its equilibrium state. Now we consider 
two superconductors, A and B, separated by an insulating barrier. If the barrier is thick enough so that 
the superconductors are isolated from each other, the time-dependent Ginzburg-Landau equation for 
each side is 



ihdt^A = H a i/j a + ki/j b , 



(2.1) 
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(2.3) 
(2.4) 



m AjB is the mass of each cooper pair, a AtB (T) and b AjB (T) are two parameters that depend on the 
temperature, k is the coupling constant for the wave-functions across the barrier. It is worth noting that 
a A>B (T) and b A)B (T) are phenomenological parameters that can be fixed in experiments. Without the 
coupling k, no mass term is generated in our paper. We will set m A = m B in the following. Actually, 
equations (12.11 ) and ( 12.21 ) are used to describe the Josephson effect when a voltage V is applied between 
the two superconductors and then one can replace the Hamiltonian with H A = eV and H B = —eV. We 
use the coupled time-dependent Ginzburg-Landau equations in order to obtain the effective acoustic 
metric. Note that the similar computation was done for a two-species Bose-Einstein condensate by 
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using two-component time-dependent Gross-Pitaevskii equations 0. One may regard 6^ and bs as 
self-interactions U aa an d Ubb m GO- But different from the two component Bose-Einstein condensate, 
here the two superconductors A and B are separated by a thin film and the interactions between the two 
superconductors are very weak. So the interactions bAB an d b B A can be neglected. After obtaining the 
metric, we will ask what we can learn from the Josephson effect for our understanding on the acoustic 
black hole physics. 

The two macroscopic wave function can be written in the form 

^A = JP~Ae ie \ (2.5) 
*!>b = VPB~e i6B . (2.6) 

The equations of motion then becomes 

d t pi H V ■ (piVi) + -JpWj sin(0i - 8j) = 0, (2.7) 

rrii h 

h 2 V 2 



where vi = — and i, j = A, B (i ^ j). The first term in the right hand of (12.81 ) corresponds 
to the quantum potential. The above two equations are completely equivalent to the hydrodynamic 
equations for irrotational and inviscid fluid apart from the quantum potential and the k term. In the long- 
wavelength approximation, the contribution coming from the linearization of the quantum potential can 
be neglected. Since the current densities p i in the superconductor A is not much different from that in 
the superconductor B, the quantum potential can be neglected in our derivation of the acoustic metric. 
The physics of the Josephson current is mainly described by the k term in equations (12.71 ) and (12.81) |fT9l . 
The Josephson relations for the pair density can be obtained from (12.71) 

/j 2k 

d t (pA - Pb) H V • (paVa - Pb$b) + —^PaPb sm(6 A - Ob) = 0. (2.9) 

V(Ia h 

It can be reduced to the standard Josephson relation when the second term is dropped out, that is to say, 

j = 2ed t (p A - Pb) = ~—\/PaPb sin(0 B - A ) =j c sin9, (2.10) 

n 

where 

4eK 



jc = —VPA~P~B~, 9 = 9 B -9 A . (2.11) 
n 

Actually, only in the Josephson junctions, the speed of electrons cannot be neglected. The kinetic 
energy of one cooper pair tunneled from one side to the other side changes little. That is why one can 
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replace the Hamiltonian with Ha = eV and Hb = —eV and ignore the kinetic terms of the Hamiltonian 
in (12.81) . Following this procedure, we obtain the second Josephson equation 

d(e B -e A) = 2eV 

dt ft v 

In the insulating barrier, the density of the current carriers pj is much lower than the density of current 
carriers in the superconductors p s . For those insulating films who share the same cross-sectional area 
with the superconductors on each side, they share the same current. That is to say 

p j v j eS 1 = p s v s eS 2 , (2.13) 

where Si = S 2 = S is the cross-sectional area, Vj denotes the speed of the current carriers in the film 
(junction), and v s the speed in the superconductors. Therefore, the speed of the current carriers in the 
film (junction) must be much bigger than their speed in the bulk: 

Vj->v s . (2.14) 

On the other hand, if the two superconductors, A and B, are not separated by an insulating barrier, but 
by the same superconducting material with a cross-sectional area much smaller than the cross-sectional 
area of each side (i.e. Si <C S 2 ), this is another kind of Josephson junctions. In this case, the density of 
current is same everywhere in the around circle (i.e. pj = p s ). Thus, we still have Vj ^> v s . These can 
justify why we can drop the second term in (12.91) and kinetic energy terms in (12.81) when we derive the 
Josephson relations and how it is possible for the creation of an acoustic black hole in the Josephson 
junction because the speed of the current carriers can be so fast that it may exceed the local speed 
of sound. In this paper, we will construct an acoustic black hole from the supercurrent tunneling by 
considering the linearized perturbations of equations (12.71) and (12.81) . 



3 The acoustic black hole metric 

Now let us consider a fixed background (p i0 , #j ) with small perturbations pi = pio+pnand9 = 9 i0 +9ii 
(i — A, B ). The leading order equations for (p i0 , 6 i0 ) can be written as 

h k 

dtp Ao H V • (paoVao) + tVPaoPbo shi(6Uo - 6bo) = 0, (3.1) 

rriA ft 

ftd t 9 A o = ~-^-Vaq ~ a A(T) - b A (T)p A o - kJ — cos(6 B0 - 6 A0 ). (3.2) 
* V Pao 

dtPBo + —V • (pboVbo) + tVPboPao sm(9 B o - Oao) = 0, (3.3) 

VTiB ft 

ftd t 8 B o = --tt^Io - a B{T) - b B (T)p B o - K \ — cos(6> B0 - Oao)- (3-4) 

& V PBO 
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Linearizing the equations (12.11) and (12.21) . we obtain the two coupled equations for the perturbation 
of the phases 



do 
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dt 

dO B i 
dt 
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(3.5) 
(3.6) 



The coupled equations for the density perturbations are given by 

fa 2k 

dtPAi + V • ( Pa<N6ai + PaiVao) = — \fp~AoPmiO bi - Oai), 

VTLa h 

fa 2k 

dtPBi + V ■ ( psoV^ + PbiVbo) = —\/PaqPbo{Qai -0 B i)- 

TUb h 



(3.7) 
(3.8) 



The above equations governing the perturbation of phases and density are very hard to decouple. We 
need impose some constraints on the background parameters for the purpose of deriving the acoustic 
metric. The superconductors A and B can be the same so that pao = Pbo = Po, ^a = &b = b, 
rriA = rriB and the background phases also can be set to be equal (6Uo = #bo)- This implies that without 
perturbations there are no currents crossing the junction and background velocities vao = v B o = Vq. But 
when the phase 9 is fluctuated, the supercurrent tunneling happens in the junction and the background 
velocity v can be regarded as the function of the space variables Xj. In the region where vq exceeds the 
"sound velocity" c s , an acoustic black hole forms. In the following, we will see how this can happen. 
The coupled equations for phase and density perturbation can be written as 

b, , k pai 

h po 



0, 



tWAl 



0bi)+v o -V(9 A i-9 



Blj 



iP 



Al 



dt{ P 



Al 



PB\, 



m 



-p V(fl 



Al 



e 



B\, 



(p 



Al 



pBl)v 



PBI, 



Pbi 



h 



Ml 



TBI 



(3.9) 
(3.10) 



It is convenient to introduce the notations 



9\ = Oai — Obi, Pi = Pay — Pbi, 



and 



X 



bpp - K 
hp 



After combining (13.91) and (13.101) as a single equation, we have the wave equation for Oi 



-d t 
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(d t 1 + v o -V0 1 ) 
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(3.11) 
(3.12) 

(3.13) 
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The above equation is comparable with a massive Klein-Gordon equation in curved space-time 
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We can therefore read off the inverse acoustic metric 

-1 



-MV^gg^d^) - m% = o 



(3.14) 
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where the local speed of sound is defined as 

„2 



hpo _ bp - k 
m m 



Note that in absence of the coupling constant k, the local speed of sound has the form 

h 



where £(T) is the Ginzburg-Landau coherence length ^(T) 
By inverting (13.15k we determine the metric 



(3.16) 



(3.17) 



(3.18) 



\/2m|a(T)|' 
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(3.19) 



In the presence of an external magnetic field, we will show in the following that the structure of the 
acoustic black hole may have "draining bathtub" form. The general acoustic metric is given by 

hp 



me* 



- (cf - ifydt 2 - 2v ■ drdt + df-df 



(3.20) 



where the corresponding horizon locates at c s = vq. 

In the cylindrical coordinate (r,9,z), suppose that the superconducting current is along the z- 
direction and the magnetic field A is along the ^-direction. In this case, the background fluid flow 
will be bended by the magnetic field. It is convenient to set v z (r, 9, z) ^ and vg(r, 9, z) ^ and let 
v r — 0. The metric (13.201) becomes 



ds z 



mc. 



Vrac,/ 



— (c 2 — v^)dt 2 — 2v* Q ■ drdt + dr ■ dr 
c 2 s dt 2 + (v g dt - rd9) 2 + (v z dt - dz) 2 + dr 2 



(3.21) 
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where v% 



v\ + Vg. If we make the coordinate transformations 



dt = dr 




(3.22) 



d9 = dd 



v d Vz 



dz, 



(3.23) 



r(c 2 s - v I) 



then the line-elements of the metric can be written as 



dg2 = rnpo\ | _ [c 2 _ {v 2 + v 2 )]dr 2 _ 2rVedrM + _^ 




(3.24) 



The formation of an acoustic black hole requires ^ > 1 in some regions. In ifTOl . the authors 
pointed out that it would be very difficult to realize acoustic black holes by using type I and type II 
superconductors. Especially, for type I superconductors, when the speed of superconducting electrons 
is equivalent to the "sound velocity", say vq = c s , the superconducting phase is broken and return to 
the normal state. It was expected to form an acoustic black hole in the region of the spiral vortex core 
in a type II superconductor. The calculation shows that £ < r < v^2£, electron velocity may exceed the 
sound velocity [[Toll , where r denotes the distance from the vortex core and £ the coherence length. But 
experimental verification of superconducting electron speed in the bulk of a superconductor is not easy. 
From (12.111) (i.e. j c = 2ev l p ^ 4enp /h) and (13.171) . we know that ^ > 1 means 



As an example, let us consider the material PbCu at T c = 7.2K with coherence length £ = 80rim |fT9l . 
Then (13.251) requires the coupling constant k > 1.788 x 10~ 18 J , which is possible in experiments. 
The Josephson junctions and related instruments, such as superconducting quantum interference device 
(SQUID), may open a door to build an acoustic black hole directly. 

4 A microscopic picture 

From (13.171) . we know that the local speed of sound in the junction could be smaller than that in the 
bulk of the superconductor for a positive-valued coupling constant k. In the above derivation, we 
have used the Ginzburg-Landau theory, which is very useful in describing qualitative and macroscopic 
behaviors. In order to have a clear picture for the formation of acoustic black holes, now we present a 
microscopic description by using the BCS theory. In fact, Gorkov in 1959 proved that the Ginzburg- 
Landau theory can be derived from full the BCS theory in a suitable limit IfTTll . The microscopic model 




(3.25) 
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can be constructed by using the tunneling Hamiltonian 



H = 



Hr + H L + H T , 

= J2( T ><p c l c p° +h - c -y 

kpcr 



(4.1) 
(4.2) 



The Hamiltonian is identical with equations (12.11 ) and (12.21) . H R is the Hamiltonian for particles on the 
right side of the tunneling junction. Similarly, Hl has all the physics for particles on the left side of the 
junction. The tunneling is caused by the term H T and T pcr denotes the tunneling matrix that can transfer 
particles through an insulating junction. The derivation of the tunneling current led to two terms: the 
single-particle terms and the Josephson term. We only consider the Josephson term here. 



ihdtipi = Hiipi + Kipj, 
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where e f is the Fermi energy. In this sense, the formation of an acoustic black hole requires 

> 1. 
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As pointed out in section 2, in the insulating barrier, the density of the current carriers is much lower 
than that in the superconductors, but the speed of the current carriers in the film should be much bigger 
than their speed in the bulk. Therefore, there may exist a region in where v < c s that continuously 
connected to the region v > c s . 



5 An acoustic black hole in the Josephson junction 

Let us consider a weak link tunnel junction with a magnetic field B x (y)i applied along the x-direction, 
as shown in Fig.l. The junction is of thickness 2a normal to the z-axis with cross-sectional dimensions 
d and w along y and x, respectively. We assume that the external magnetic field is larger than the 
field produced by the currents. The applied field is derived from the vector potential A = B x (y)k. In 
the barrier film the material is normal and the magnetic field is constant- valued B x (y) = B , but the 
magnetic field decays exponentially into the superconductors on either side of the junction. 
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From (13.211) . we know that the magnetic field can change the direction of the fluid flow. The for- 
mation of an acoustic black hole should satisfy the condition v' l > c s . Note that when we derive the 
acoustic metric, we consider the perturbations around a fixed background (pio, Ow) without the fluctua- 
tions of the magnetic field. The vector potential A is regarded as an external source. We know that the 
phase 9q is determined by the magnetic field and the supercurrent. Let us first review the derivation of 
the Josephson junction diffraction equation^ see lfT8l[T9l for more details). Consider a rectangle circle 




Figure 1: Application of a magnetic field B x transverse to the Josephson junction. An acoustic black 
hole may be created in the junction. 

PC1QC2P in the Josephson junction, where P and Q locate at the middle of the junction. We can 
neglect the thickness of the film and assume Ay deep inside the superconductor where the induced 
current decays away (i.e. the magnetic field is vanishing). The integration of \ z ■ dl along the rectangle 
circle PCiQC 2 P is zero. Therefore, only the magnetic field contributes to the change of the phase 

W = ^A. (5.1) 

n 

In the superconductor A, the integral around C\ path gives 

OQl (y) - 9 0Pl (y + Ay) = ^ / A- dl. (5.2) 

For the C 2 path 

0op 2 (y) - OQ2 (y + Ay) = ^ / A ■ dl. (5.3) 

h JC2 
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Then we have 

9 (y + Ay) - 6 (y) = ^^A-dl. (5.4) 
By using the Stokes theorem, we find that 

= ^ J J B-ds=jB x (2K + 2a)Ay. (5.5) 

Now, we have 

2ck 

V y 6 = —B x (y), (5.6) 

h 

where A = 2k l + 2a denotes the effective thickness of the junction and k l is the penetration depth of 
the magnetic field. The phase 6 then depends on the coordinate y 

9 (y)^ 2 -^y + c 1 . (5.7) 

It is worth noting that we can replace B x with the total flux across the junction $ j = AdB x , 

2eAB x 2ttMB x 2tt$j 
fd W dfo ' 



where O = 2~ is the magnetic quantum flux. If this is substituted in Eq. (12.10l) and integrated over the 
area S = wd of the junction, we have 

sin(^) 

Is(B)=j c S Jl sin Cl . (5.9) 

00 

We call this the Josephson junction diffraction equation. This equation indicates that the nth maximum 
of the current I s occurs at the flux value $j = (n + \)4>o, but cancels for $j = n0 o , where n is an 
integer. 

We are interested in (15.61) . since this equation gives us v y = ^V y 0o = ^B x (y). Therefore, in the 
cartesian coordinate system, from (13.191) we have the metric 



ds 2 = (— — ^ — (c 2 . — (v 2 + Vy))dt 2 — 2v z dzdt — 2v y dydt + dxi ■ dxi 
Taking the coordinate transformation 



(5.10) 



dt = dr -dy V — -dz, (5.11) 

cj-v 2 - v 2 c 2 s -v 2 - v 2 



we obtain 

2 f hp n r / 



ds 



) - \ c l - ( v z + vl)]dr 2 + [ 5ij + — — ^ —dx l dx j 

\ ci - (v 2 + v 2 ) 



(5.12) 
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where we should note that c s 



bpp—K 



¥ and v y 



—B x {y). Note that the transformation 
from the metric (15.101) to (15.121) means that the resulting metric (15.101) is a static one. The variable t 
measures the time of the background fluid and r is a redefined time. It would be non-trivial to consider 
the case in which the total phase change across the junction is 2mr, with n Josephson vortices side 
by side in the junction, each containing one flux quantum. Then the total current in the junction is 
vanishing because the supercurrent flows down across the junction on the left and up on the right. The 
current flows horizontally within a penetration depth A inside the superconductor to form a closed loops. 
These current loops encircle flux and the resulting configuration is known as a Josephson vortex. The 
supercurrent along the z-direction can somehow be regarded as a constant value and then v z does not 
change in the junction. 

As a special condition, let us consider the Hawking temperature contributed by v y . In the vicinity of 
the horizon, we can split up the fluid flow into normal and tangential components (i.e. v = v 2 + \ y ) and 
we choose j as the unit vector field that at the horizon is perpendicular to it |[2). From the definition, we 
know that the horizon locates at c s — v y . In this case, the Hawking temperature at the event horizon is 
given by 



2nk f 



d y (c s 



horizon 



(5.13) 



More explicitly, 



T, 
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(1.2 x 10~ 9 fTm) 
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d y (c s 



DyCg 
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2eAfi j z 
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horizon 



(5.14) 



where we have used the Maxwell equation V x B = /i J. If the speed of sound is to be a position- 
independent constant, the resulting Hawking temperature then has the form 



T H = (1.2 x WT*Km) ( 



( 1 ) 


2eA/j j z 


\ 1000ms- 1 J 


m 



horizon 



(5.15) 



This is a very interesting result which indicates that for big enough value of the effective thickness of 
the junction A and the tunneling current j z , the Hawking temperature would be detectable in the future. 
As an example, let us estimate the Hawking temperature of a particular kind of Josephson junction: 
given the tunneling current j z = 5 x 10 7 A/m 2 and the effective length A = 50nm (including the 
penetration depth) (TH), the resulting Hawking temperature is about Th ~ 10~ 7 Kelvin. This value 
varies for different tunneling currents j z . Compared with the Hawking temperature (~ nano Kelvin) 
of acoustic black holes in Bose-Einstein condensate, the temperature obtained here is two orders of 
magnitude higher that maybe possible for the future experiments. 
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6 Conclusion 



In summary, we have presented a version of acoustic black holes by using the Josephson effect. We 
started from two coupled Ginzburg-Landau equations with a coupling constant k for Josephson junc- 
tions and reviewed the basic equations for the Josephson effect. 

The acoustic black hole metric can be obtained from the perturbation equation for 9 1 and p\. The 
advantage of creating acoustic metric by using the Josephson effect is that the coupling constant k can 
be tuned. So that the sound velocity- c s can be tuned to be very small and then v 0i > c s would become 
easier. We discuss in detail how to describe an acoustic black hole in the Josephson junction and 
write the metric in the language of the superconducting electronics. Finally, we estimate the Hawking 
temperature of acoustic black hole created in the Josephson junction. Although we have set up a 
theoretical model for acoustic black holes in Josephson junctions, the experimental detection of the 
Hawking temperature would be difficult. Our result indicates that the Hawking temperature strongly 
depends on the tunneling current. The enhancement of the Josephson current is thus crucial for the 
measurement of the acoustic black hole. 

On the other hand, the obstacle of detecting Hawking radiation may come from its instability against 
other mechanisms. For instance, for type I superconductors, when v = c s , the superconducting phase 
is broken and return to normal states ifTOll . In EOl . the authors studied acoustic horizons in the quantum 
de Laval nozzle. They solved the Gross -Pi taevskii equation and found that both in hydrodynamic and 
non-hydrodynamic regimes there exist dynamically unstable regions associated with the creation of 
positive and negative energy quasiparticle pairs in analogy with the gravitational Hawking effect. In 
this paper, we may suffers the same problems since the Ginzburg-Landau equation and the Gross- 
Pitaevskii equation are very similar. Also, the quasinormal modes analysis of the obtained acoustic 
black holes may reveal that in the high momentum regime the configuration would be unstable against 
perturbations. We leave discussion of the quantum instability to a future publication. 
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